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Abstract 

A  novel  approach  for  solving  the  DEA  linear  programming  problems  using  a  primal- 
dual  interior-point  method  is  presented.  The  solution  found  by  this  method  satisfies 
the  Strong  Complementarity  Slackness  Condition  (SCSC)  and  maximizes  the  product 
of  the  positive  components  among  all  SCSC  solutions.  The  first  property  is  critical  in 
the  use  of  DEA  and  the  second  one  contributes  significantly  to  the  reliability  of  the 
solution. 
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1  Introduction 


Data  Envelopment  Analysis  (DEA)  is  a  linear  programming  application  from  the  economic 
and  management  sciences  where  a  desired  solution  of  the  underlying  linear  programming 
model  is  one  satisfying  the  strong  complementarity  slackness  condition  (a  SCSC  solution). 
The  set  of  SCSC  solutions  can  be  characterized  as  the  set  of  solutions  where  the  number  of 
nonzero  variables  is  maximal,  therefore  SCSC  solutions  are  not  vertex  (or  basic)  solutions 
unless  the  solutions  for  both  the  primal  and  dual  linear  programming  problems  are  unique. 
Thus,  the  standard  simplex-type  algorithms  are  generally  not  appropriate  for  computing 
SCSC  solutions. 

In  this  work  we  propose  a  primal-dual  interior-point  algorithm  for  solving  the  DEA  linear 
programming  problems.  The  solution  found  by  this  interior-point  algorithm  is  not  only  a 
SCSC  solution  but  it  is  also  the  center-most  solution  in  the  sense  that  it  maximizes  the 
product  of  the  positive  components  among  all  the  SCSC  solutions  (the  so-called  analytic 
center  of  the  solution  set).  These  properties  are  crucial  in  the  use  and  reliability  of  DEA 
since  among  those  SCSC  solutions  it  is  desirable  to  compute  one  that  is  in  some  sense  as  far 
away  as  possible  from  the  relative  boundary,  avoiding  in  this  way  that  small  changes  in  the 
data  may  affect  the  DEA  results  (see  Thompson,  Dharmapala,  and  Thrall  [16]  ). 

This  paper  is  organized  as  follows.  In  the  next  section  we  introduce  DEA  and  give  a 
formal  definition  of  the  SCSC  solutions.  Section  3  contains  fundamental  background  for  the 
proposed  algorithm.  In  Section  4  we  present  the  algorithm  and  we  study  its  behavior  in 
Section  5.  Numerical  results  considering  real  world  DEA  problems  are  included  in  Section 
6.  Finally,  in  Section  7  we  give  some  concluding  remarks. 
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2  Data  Envelopment  Analysis  (DEA) 


DEA,  as  presented  in  Charnes  and  Cooper  [2],  is  a  technique  for  measuring  the  amount 
of  relative  efficiencies  and  inefficiencies  in  a  production  system  characterized  by  data  units 
called  Decision  Making  Units  ( DMU ). 

Each  Decision  Making  Unit  DMU  is  a  vector 

Rm,(y,  x)  >  0  and  neither  y  nor  x  is  zero. 

The  amounts  of  output  r  and  input  i  for  DMUj  are  represented  respectively  by  yrJ  and 


y 

—x 


£  R,+m  where  y  £  R*,x  € 


x 


•r 


A  DEA  Data  Domain  D  considers  a  set  of  n  decision  making  units  DMU\ , . . .  ,DMUn 
and  a  collection  of  technologies  which  is  characterized  by  a  data  matrix  P  —  ( P\...Pn )  = 

y  \ 

£  R(m+s)xn  with  p.  =  DMUj,  Y  =  (j/i-.j/n)  £  RsX"  and  a:  =  (*!...*„)  £  Rmxn. 

~x }  .  , 

A  standard  assumption  for  DEA  is  that  no  two  columns  of  P  are  proportional. 

The  production  possibility  set  associated  to  the  data  matrix  P  is  defined  as 


lC(V)=\p  = 


y 

—  X 


£  Rs+m  :  p  <  P A  for  some  A  >  0,  A  £  Rn,  and  (y,  x)  >  0  >  . 


Let  DMU0  = 


Vo 


denote  one  of  the  DMU' s  (i.e.  o  £  {1, . . . ,  n}  ). 


\  / 

Consider  the  linear  programming  problem  associated  to  DMU0 


minimize 
subject  to 


e 

S  =  P\- 
A  >  0, 


Vo 

-0xo 


>0, 


(1) 


where  9  £  R,  A  £  R"  and  S  £  Rs+m. 

Problem  (1)  is  referred  to  as  the  envelopment  (or  primal)  program  for  DMU0- 

This  problem  is  the  CCR  ratio  model  introduced  by  Charnes,  Cooper,  and  Rhodes  (1978). 
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The  present  paper  focuses  attention  on  this  DEA  model  although  our  analysis  can  be  ex¬ 
tended  to  the  BCC  linear  model  (Banker,  Charnes,  and  Cooper  (1984)). 

Let  (Q*,\l,S*)  be  a  solution  of  (1).  The  DEA-radial-efficiency  for  DMU0  is  defined  as 
the  optimal  value  6*. 

The  DMU0  is  said  to  be  DEA-radial-efficient  if  6*  —  1  and  the  DMUa  is  said  to  be 
DEA-radial-inefficient  if  0*  <  1  .  This  definition  is  equivalent  to  saying  that  a  DMU  is 
DEA-radial-efficient  if  it  lies  in  the  frontier  of  the  production  possibility  set. 

The  definition  of  efficiency  used  in  DEA  is  based  on  the  idea  that  to  remain  in  the  pro¬ 
duction  possibility  set  a  DMU  is  technical  efficient  if  the  inputs  and  outputs,  corresponding 
to  that  DMU ,  can  not  be  respectively  decreased  nor  increased.  Therefore,  a  DMU0  can  be 
a  boundary  point  ( 6 *  =  1)  and  not  be  efficient.  For  this  reason,  a  given  DMU0  is  called 
DEA-efficient  if  0*  =  1  and  S*  =  0  for  all  the  solutions  (0*,  S*,  A*)  of  problem  (1). 

There  is  an  alternative  dual  approach  to  DEA-radial-efficiency  for  DMU0  using  the  mul¬ 
tiplier  space  W0  defined  as 

W0  =  {w  =  (  'j  G  Rs+7n  :  ha(w )  >  hj(w)  for  all  j  =  1, . . .  ,n,  w  >  0  and  u^0,r/0}. 


Here  the  functions  hj 


are  output/input  ratios  defined  as 


hj(w)  =  uTyj/vTxj  for  j  -  1, . . . ,  n , 


u 


for  each  w  =  ,  u  G  Rs,  v  G  Rm,  w  >  0  such  that  vTxj  >  0. 


Let 


W”  =  {u;  G  W0  :  vTx0  =  uTya  =  1}. 


The  set  WTan  is  called  the  normalized  multiplier  set  for  DMU0  ■ 
Observe  that  W”  ^  0  if  and  only  if  W0  7^  0. 

Consider  the  linear  programming  problem 
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maximize  uTy0 
subject  to  vTxa  =  1, 

-t,  =  uTy i  -  vTx j  =  wTPj  <  0,  j  =  1 , . . . ,  n  (2) 

where  u  €  Rs,  v  G  Rm  and  t  €  Rn- 

Problem  (2)  is  the  dual  problem  of  (1)  and  it  is  referred  to  as  the  multiplier  (or  dual) 
program  for  DMUa. 

It  is  easy  to  see  that  (w*,t*0)  is  a  solution  of  problem  (2)  if  and  only  if  w*  <E  W™.  Therefore, 
by  duality  theory,  W0  7^  0  if  and  only  if  9*  =  1.  Hence,  DMU0  is  DEA-radial-efficient  if 
W0  0  and  is  DEA-radial-  inefficient  if  W0  =  0. 

Let  be  the  set  of  all  the  DEA-radial-efficient  DMUj  and  J\f  be  the  set  of  all  the 

DEA-radial-inefficients  DMUj ,  for  j  =  1, . . .  ,n. 

The  DEA-radial-efficiency  can  be  computed  using  either  of  the  two  linear  programs  (1) 
and  (2).  Hence,  the  solutions  of  these  problems  answer  the  question  of  DEA-radial-efficiency. 
However,  among  those  DMU's  that  are  DEA-radial-efficient  (or  DEA-radial-inefficient), 
there  are  important  differences  depending  on  their  multiplier  sets.  Thus,  the  set  of  DMU  s 
may  be  partitioned  into  six  classes  (see  [3]  and  [4]  for  more  detail): 

E  =  {DMUj  e  7 Z£  :  dim  Wj  =  s  +  m}, 

E'  =  { DMUj  e  7 ZS  :  dim  W j  <  s  +  m  and  exist  u;  >  0, u>  €  Wj}, 

F  =  {DMUj  €  7l£  :  every  w  €  Wj  has  at  least  one  zero  component), 

NE  =  {DMUj  e  AT  :  DMU]  €  E), 

NE'  =  {DMU,  6  N  :  DMU ]  €  E'}, 

NF  =  { DMU j  €  AT  :  DMU ]  €  F}. 
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Here  DMU'}  =  P3(0*)  = 

\  ~@jxj  ) 

The  elements  of  E ,  E',  and  F  are  called,  respectively,  DEA-extreme  efficient,  DEA-non- 
extreme  efficient,  and  DEA-weak  efficient.  It  can  be  seen  [4]  that  any  given  DMU0  is 
DEA-efficient  if  it  belongs  to  E  U  E' ,  or  equivalently,  if  there  exists  a  positive  multiplier 
vector  w  £  W0. 

The  classification  and  characterization  of  the  DMU's  play  a  very  important  role  in  DEA. 
Charnes,  Cooper,  and  Thrall  [4]  present  a  structure  for  this  characterization  using  the  linear 
programming  problems  (1)  and  (2).  A  very  interesting  result  from  [4]  is  Lemma  10A  which 
gives  a  way  to  classify  the  different  DMU’s,  but  requires  the  computation  of  a  solution 
{0oAo,So,wo,to)  of  problems  (1)  and  (2)  satisfying 

SIWo  =  fro  =  0 

and 

S0  +  w0>  0,  t0  +  A0  >  0. 

These  solutions  are  said  to  satisfy  the  strong  complementarity  slackness  condition  (or  are 
SCSG  solutions).  It  is  well  known  (see  Spivey  and  Thrall[14]  for  a  proof)  that  there  always 
exists  a  SCSC  solution  for  any  linear  programming  problem  which  is  both  primal  and  dual 
feasible. 

I  he  power  of  the  SCSC  approach  depends  on  how  far  wa  is  from  the  boundary  of  the 
set  W0.  In  order  to  explain  this  idea  (see  Thompson,  Dharmapala,  and  Thrall  [16]  for  more 
details)  let  us  consider  a  DMU0  being  DEA-extreme  efficient  (i.e.,  DMUa  £  E ).  Then,  for 
every  SCSC  solution  (0O,  Xa,  Sa,  wa,  tQ),  w0  is  the  center  of  an  (m  +  s)-dimensional  sphere 
(hypersphere)  which  is  a  subset  of  W0  and  for  every  point  w'a  in  this  sphere,  h0(iu'0 )  >  hj(w'0 ) 
for  all  j  ^  o.  Therefore,  one  has  the  reasonable  expectation  that  DMU0  remain  DEA- 
extreme  efficient  for  relative  small  changes  in  the  data.  Moreover,  if  we  define  the  sensitivity 
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function 


d0(wo)  =  h0(w0)  -  hk(w0),  with  hk(w0)  =  n?a xhj(w0) 

we  expect  the  DEA  results  to  be  less  sensitive  to  changes  in  the  data  when  d0(w0)  is  larger. 
Therefore,  for  DEA  purposes,  the  best  choice  among  all  the  SCSC  solutions  is  the  one  that 
maximizes  the  function  dQ(w )  over  the  set  of  dual  multipliers  w  with  (0,  A,  S,w,t)  a  SCSC 
solution.  Let  us  call  this  solution  the  DEA-center  solution  ( 0co ,  Aqj  S £,  wc0,  and  dca  =  d0{wc0). 

The  DEA-center  is  an  useful  theoretical  definition  for  understanding  the  idea  of  centrality 
that  best  fits  the  DEA  requirements,  however  its  practical  use  depends  on  the  development 
of  an  effective  numerical  method  for  computing  it.  In  this  paper  we  propose  an  algorithm 
that  computes  a  SCSC  solution  called  the  analytic  center.  This  solution  satisfies  a  notion  of 
centrality  that  we  hope  is  close  to  the  idea  of  centrality  given  by  the  DEA-center,  however 
the  exact  relationship  between  the  analytic  center  and  the  DEA-center,  as  well  as  the  com¬ 
putation  of  the  DEA-center  solution,  are  important  topics  that  will  not  be  covered  by  this 
paper  and  deserve  further  study. 

3  Algorithm  Background 

Simplex-type  methods  consider  the  vertices  of  the  feasibility  set  of  the  linear  programming 
problem  and  the  final  solution  is  an  optimal  vertex.  On  the  other  hand,  interior-point 
methods  generate  an  infinite  sequence  of  points  in  the  interior  of  the  feasibility  set  and  the 
algorithm  stops  when  the  iterates  are  sufficiently  close  to  a  solution. 

There  are  several  different  classes  of  interior-point  methods  for  linear  programming  (see 
Gonzaga[7]  for  a  survey).  In  this  work  we  are  interested  in  primal-dual  interior-point  al¬ 
gorithms.  We  introduce  in  this  section  some  theoretical  background  concerning  these  algo¬ 
rithms  as  well  as  a  discussion  of  the  main  results  necessary  to  motivate  the  development  of 
our  algorithm. 
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We  employ  the  standard  notation  used  in  the  literature  concerning  primal-dual  interior- 
point  algorithms. 

Consider  the  linear  programming  problem  in  the  standard  form: 


minimize  cTx 

(3) 

subject  to  Ax  =  b,  x  >  0, 

where  c,  x  €  Rn,  b  G  Rm,  A  €  RmXn(m  <  n)  and  A  has  full  rank  m.  Associated  with 
problem  (3)  is  the  dual  linear  program 

maximize 
subject  to 

where  y  G  Rm  and  zGR". 

Let 

?  =  {(z,  V,  z):  Ax  =  b,  AT y  +  z  =  c,  (x,  z)  >  0}. 

A  point  (x,y,  z)  G  is  said  to  be  strictly  feasible  if  x  and  z  are  strictly  positive. 

It  is  well  known  that  x  is  a  solution  for  (3)  and  ( y ,  z )  is  a  solution  for  (4)  if  and  only 
if  (x,y,z)  G  T  and  the  duality  gap  cTx  —  bTy  =  0,  or  equivalently,  xTz  =  0.  This  latter 
condition  is  called  the  complementarity  slackness  condition. 

Primal-dual  interior-point  algorithms  attempt  to  solve  the  primal  and  dual  linear  pro¬ 
grams  simultaneously  by  generating  a  sequence  of  strictly  feasible  points  ( Xf.,yk,Zk )  such 
that  the  duality  gap  xj^z*,  converges  to  zero. 

Observe  that  xT z  =  0  with  (x,z)  >  0  if  and  only  if  x,z,  =  0  for  all  i  =  1, . . .  ,n.  The 
optimality  conditions  for  (3)  and  (4)  can  be  written  as 


bry 

ATy  +  z  =  c,  z  >  0, 


F(x,y,z) 


( 


\ 


Ax  —  b 
ATy  +  z  -  c 
XZe 


=  0,  (x,z)>0, 


(5) 
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where  X  =  diag(x),  Z  =  diag(z)  and  e  =  (1, 1, . . . ,  1,  l)r  G  Rn. 

Most  primal-dual  interior-point  algorithms  are  variants  of  Newton’s  method  applied  to 
problem  (5). 

We  denote  the  solution  set  of  problem  (5)  by 


S  =  {(x,y,z)  :  F(x,y,z)  =  0,(x,z)  >  0}. 


A  point  ( x,y,z )  in  S  is  said  to  be  a  strict  complementarity  solution  if  x,-  +  >  0  for  all 

i  =  1, . . . ,  n. 

For  a  given  DMUa ,  the  DEA  problems  (1)  and  (2)  discussed  in  the  previous  section  can 
be  written  in  the  standard  forms  (3)  and  (4)  by  defining 


A  = 


-In 


-/.XI 


0  P 

Xn 


\ 


€  r(™+*)*(™+"+*+i))C  _  ]  g  R 


m+n+s+l  1  _ 


,b  = 


GR 


m-fs 


(S\ 

!  w\ 

Therefore,  x  = 

6 

G  Rm+n+s+1 ,  and  z  — 

0 

\t; 

G  R,m+n+5+1 


In  this  work  we  assume  the  following  standard  conditions  for  (3)  and  (4): 


(Al)  {x  G  R"  :  Ax  =  b,  x  >  0}  ^  0,  and 
(A2)  {(y,  z)  G  Rm+n  :  ATy  +  z  =  c,  z  >  o}  ±  0. 

We  are  particularly  concerned  with  the  case  when  5  is  not  a  singleton  set,  i.e., 


ri(S )  7^  0 

where  ri(S )  denotes  the  relative  interior  of  S.  In  this  case,  both  problems  (3)  and  (4)  do 
not  have  unique  solutions. 

Under  the  previous  assumptions,  the  solution  set  S  has  the  following  interesting  structure: 
(i)  S  7^  0  is  bounded  ;  (ii)  all  points  in  the  relative  interior  are  SCSC  solutions  (and  all  points 
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on  the  relative  boundary  are  not);  (iii)  the  zero- nonzero  pattern  of  points  in  the  relative 
interior  is  invariant.  See  Zhang  and  Tapia  [19],  Charnes  et  al.  [4]  and  El-Bakry  et  al.  [5]  for 
proofs. 

Therefore,  for  any  (x*,y*,z*)  £  n(5),  the  following  index  sets 

It  =  { i  :  x*  >  0, 1  <  i  <  n}  and  J+  =  {i  :  z\  >  0, 1  <  i  <  n} 
are  independent  of  the  choice  of  (x*,y*,z*).  Moreover,  by  strict  complementarity 

4+U4+  =  {l,2,...,n}  and  I+f\I+  =  t 

Because  of  this  structure  of  the  solution  set  S  the  solutions  in  ri(S)  may  be  characterized 
as  those  for  which  the  number  of  nonzero  components  is  maximal.  Among  them,  there  is 
a  solution  that  may  be  thought  of  as  the  center  solution  in  the  sense  that  it  maximizes 
the  product  of  the  positive  components.  This  solution  is  called  the  analytic  center  of  the 
solution  set  and  it  was  introduced  by  McLinden  [11]  in  a  general  setting  and  independently 
by  Sonnevend  [13]  in  the  context  of  linear  programming.  It  is  interesting  to  mention  that 
Goldman  and  Tucker  [6]  introduced  a  concept  of  central  solution  for  game  theory  which  may 
be  related  with  the  concept  of  analytic  center  for  linear  programming.  However,  the  idea  of 
centrality  defined  by  the  analytic  center  is  not  exactly  the  same  as  the  central  solution  of  a 
game  as  defined  in  [6]. 

Formally,  the  analytic  center  of  the  solution  set  S  is  defined  as: 

(x*,y*,z*)  =  argmax{V>(:r,2)  :  (x,y,z)  €  5}  (6) 

where 

0(x,z)  =  n  xi  II  2*- 

>e/+  .€/r+ 

Equivalently,  in  (6)  one  can  replace  ip(x,z)  by  its  logarithm,  i.e., 

In  ip(x,  z)=  J2  In  Xi  +  Yt  ln  *.-• 

•€4+  «e/+ 
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The  central  path  of  problem  (5)  is  defined  as  the  set 

Vc  =  {(x(y),y(y),z(y))  is  strictly  feasible  :  X(y)Z(y)e  =  ye},  for  y  >  0. 

This  is  equivalent  to  saying  that  a  strictly  feasible  point  (x,  y,  z)  is  on  the  central  path  if  and 
only  if  it  satisfies  xxz\  —  x2z2  =  ...  =  xnzn. 

These  notions  of  analytic  center  and  central  path  are  well-defined  under  the  assumptions 
(Al),  (A2).  For  more  details  see  McLinden  [11],  Megiddo  [12]  and  Gonzaga  [7]. 

A  very  interesting  result  from  McLinden  [11]  is  Theorem  9.  In  the  case  of  linear  pro¬ 
gramming,  it  states  that  the  central  path  intersects  the  solution  set  at  the  analytic  center, 
i.e.,  the  central  path  point  (x(y),  y(y),  z(y))  converges  to  the  analytic  center  (x*,y*,z*)  as 
y  converges  to  zero.  See  also  Proposition  8.2  in  Megiddo  [12]  and  the  discussion  preced¬ 
ing  it.  This  result  has  played  an  important  role  in  the  development  of  most  primal-dual 
interior-point  algorithms  which  attempt  to  follow  the  central  path. 

The  interior-point  algorithm  proposed  in  this  work  has  been  designed  for  effectively  com¬ 
puting  the  analytic  center  of  the  solution  set  in  linear  programming. 

4  Algorithm 

The  first  primal-dual  interior-point  method  for  linear  programming  was  constructed  by  Ko- 
jima,  Mizuno  and  Yoshise  [9],  based  on  the  work  of  Megiddo  [12].  The  following  general 
framework  captures  the  essence  of  this  algorithm  and  contains  the  majority  of  the  primal- 
dual  interior-point  methods  that  can  be  found  in  the  literature. 

Algorithm  4.1  (Generic  Kojima-Mizuno- Yoshise  primal-dual  algorithm) 

Given  a  strictly  feasible  point  (x°,  y°,  z°).  For  Ar  =  0,1,2,...,  do 

(1)  Choose  <rk  £  [0, 1)  and  set  yk  =  ak(xk)T  zk /n. 
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(2)  Solve  the  following  system  for  (Axk ,  Ayk ,  Azk) : 


F'(xk,yk,zk) 

Ay 

=  —F(xk,yk,zk)  +  pk 

X""" 

O  O 

\ ) 

0) 

(3)  Choose  rk  €E  (0, 1)  and  compute  the  step-length  ak  =  min(l, Tkak),  where 

£k  _  _ 2ll _ 

mm({Xk)-1Axk,(Zk)-'Azk)' 

(4)  Form  the  new  iterate 

(xk+\yk+\zk+1)  =  (xk,  yk,  zk)  +  a*(Arfc,  Ayk ,  A zk). 

The  algorithm  proposed  in  this  work  is  a  modification  of  the  Kojima,  Mizuno  and  Yoshise 
algorithm,  and  is  similar  to  one  presented  in  Gonzaga  [7]. 

The  following  lemma  was  proved  in  Zhang  and  Tapia  [19].  It  provides  a  sufficient  con¬ 
dition  for  a  strictly  feasible  sequence  {(xk ,yk ,  zk)}  to  converge  to  the  analytic  center  of  the 
solution  set. 

Lemma  4.1  (Zhang- Tapia)  Let  {(xk,  yk,  zk)}  be  a  strictly  feasible  sequence  and  e  = 
(1, . . . ,  1)T  €  R”.  Assume 

ph  —  ( xk)Tzk /n  — ►  0 

and 

II Xkzk/pk  -  e ||  -*  0.  (7) 

Then  {(xk,yk,zk)}  converges  to  the  analytic  center  of  the  solution  set. 

Our  approach  to  constructing  an  algorithm  for  computing  the  analytic  center  is  to  try  to 
enforce  condition  (7). 

The  proposed  algorithm  will  utilize  two  key  ingredients: 
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•  The  notion  of  /^-neighborhoods  of  the  central  path,  defined  for  some  ft  6  (0, 1/4]  by 

M{P)  =  {(x,y,z)  :  {x,y,z)  €  T,  \\Xz  -  ye\\  <  fifi},  (8) 

where 

y  —  x  z/n. 

•  The  use  of  a  damped  Newton  method  for  solving  the  system 

Ax  —  b 

Fn(x,y,z)=  ATy  +  z-c  =0, 

^  XZe  —  y 

with  y  >  0,  in  order  to  get  to  the  central  path.  A  line-search  technique  is  used  as  a 
globalization  of  Newton’s  method.  For  each  y  >  0,  the  merit  function  considered  is 
defined  by  the  function 

/„  ■  Rn+m  ->  R 

with 

fn(x,y,z)  =  ||^x(i,  y,z)/y\\l. 

Observe  that 

f^(x,y,z)  =  \\(Xz-  ye)/y\\],  for  {x,y,z)  €  T. 

The  function  measures  closeness  to  the  central  path. 

The  idea  of  the  algorithm  is  to  consider  a  subsequence  of  gradually  shrinking  neighbor¬ 
hoods  of  the  central  path  with  /3kj  — >  0.  The  algorithm  generates  an  iteration 

sequence  {(xk,yk,zk)}  such  that  a  subsequence  (xk>  ,yk> ,zk>)  belongs  to  J\f(/3k- >).  This  sub¬ 
sequence  is  obtained  by  considering  a  sequence  of  y's  and  for  each  fixed  y  using  a  damped 
Newton  method  (with  a  line  search  globalization  strategy)  to  solve  the  system  F„(x,  y ,  z)  =  0. 
A  description  of  the  algorithm  follows: 
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Algorithm  4.2  (Long-Step  Shrinking-Neighborhood  Algorithm) 

Given  a  strictly  feasible  point  w°  =  (x°,y°,z0),  a0  G  (0,1)  and  f3°  G  (0,1). 
Choose  T]  G  (0, 1/2). 

Do  until  convergence 

(0)  Set  k  —  0. 

(1)  Define  parameters:  Set  p  =  a°(xk)Tzk/n  and  fi  =  /3k. 

(2)  Closeness  to  the  central  path:  If 


\\XkZke/p-e\\2<p 

go  to  (6). 

(3)  Compute  new  interate: 

(3.1)  Solve  the  following  system,  for  (A rk  A yk  Att^)  • 

F'(xk,yk,zk) 


fAIl 

V 

Ay 

=  -F(xt,yl‘,zk)+l, 

0 

< e  > 

(3.2)  Choose  rk  G  (0, 1)  and  compute  the  step-length  ak  =  min(l,rfcdfc),  where 


ak  = 


-1 


min((Afe)-1  Axk,  ( Zk)~iAzk ) ' 


(3.3)  Form  the  iterate 

Wk+1  =  (xk+\yk+1,zk+1)  =  ( xk ,  yk,  zk)  +  ak{ Arfc,  V,  A**). 

(4)  Line  search: 

(4.1)  IfUivk+1)  <  f^(wk)  -  (vak/2p2)^(wk)^{wk)  go  to  (5). 


(9) 
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(4-2)  If  not,  reduce  ak  and  form  the  iterate 

wk+i  =  (xfc+i5  z*+i)  =  (a>  **)  +  <*fc(Arfc,  V,  Azfc). 

(4-3)  Go  ajflin  to  step  (4-1). 

(5)  Set  k  =  k  +  1  ,  f3k  =  /3  and  go  to  (2). 

(6)  Set  p  =  a°(xk)T  zk /n. 

(7)  Do  (3.1),  (3.2),  (3.3). 

(8)  Decrease  neighborhood:  Choose  f3k+1  <  (3k . 

(9)  A;  =  k  +  1  and  go  to  (1). 

5  Algorithm  Behavior 

In  most  primal-dual  interior-point  methods  the  search  directions  used  are  linear  combinations 
of  two  fundamental  directions,  the  so-called  affine  scaling  or  pure  Newton  direction  (direction 
in  which  the  gap  is  decreased)  and  the  centering  direction  (direction  towards  the  central 
path).  This  is  also  the  case  in  our  algorithm. 

In  this  section  we  will  describe  how  the  proposed  algorithm  behaves  with  respect  to  the 
two  objectives  of  decreasing  the  gap  and  centering  the  iterates. 

The  following  propositions  address  these  issues.  Their  proofs  are  straightforward  (see 
Gonzalez-Lima  and  Tapia  [8]  for  details). 

Proposition  5.1  Let  {(xfc,  yk,  zk )}  and  {afc}  be  generated  by  the  Long-Step  Shrinking- Neighborhood 
Algorithm  (Algorithm  4-2)  and  consider  the  sets 

Icp  =  {k  >  0  :  (xk,yk,zk)  satisfies  condition  (9)  in  the  algorithm  },  i.e. 
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Icp  =  {k>  0  :  || XkZke/p  -  e||2  <  f3k},and 

I?p  =  {k>l  :  (xk~x,  yk~l ,  z*-1)  satisfies  condition  (9)  in  the  algorithm,  }  U  {0}. 

For  each  k  we  denote  by  kfp  the  largest  element  in  /+  such  that  kfp  <  k  and  by  kcp  the 
smallest  element  in  Icp  such  that  k  <  kcp  . 

Define 

=  "7  II  (1  -  o,"(l  -  <7™)),  for  /+ 

m=fc+, 

and 

ok  =  £T°,  for  k  e  ICp  U  /+. 

Then 

(i)  6  (0, 1],  for  all  k  —  0, 1, .. 

(ii)  ak  =  1  if  a*-1  =  1,  for  k  #  Icp(J  /+ 

.  <r*:i+1  <  crfc‘*+2  <  ...  <  crkcp  <  1. 

Proposition  5.2  Let  {(xk ,yk , zk)}  and  {ak}  be  generated  by  the  Long- Step  Shrinking- Neighborhood 
Algorithm  (Algorithm  f.2)  and  let  crk  for  all  k  >  0  6e  defined  as  in  Proposition  5.1. 

Then 

(i)  For  all  k  =  0,1, p  =  ok(xk)Tzk/n 

(ii)  (. xk+l)Tzk+ 1  =  (1  -  ak(l  -  ok))(xk)Tzk. 

Therefore, 

0 xk+l)Tzk +1  <  {xk)Tzk,  for  all  jb  =  0,1,2, 

If =  1)  then  (xk+1)Tzk+1  —  (xk)T zk  (the  gap  does  not  decrease).  If  ak  =  0,  then 
(xk+1)T zk+1  <  ( xk)T zk  ( maximum  possible  gap  decrease  in  that  iteration). 
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Proposition  5.3  Let  sk  =  (Axk,Ayk,A zk)  be  the  step  direction  computed  by  the  Long-Step 
Shrinking- Neighborhood  Algorithm  (Algorithm  f.2),  with  P  >  0  fixed,  and 

ffi(x,y,z)  =  \\(Xz  -  pe)/ p\\\. 

Then 

sk  is  a  descent  direction  for  f^(x,y,z)  at  (xk,yk,zk). 

The  previous  propositions  state  that  while  condition  (9)  (closeness  to  the  central  path) 
in  Algorithm  4.2  is  not  satisfied,  p  will  be  left  fixed  and  the  dominant  direction  will  be  the 
centering  direction  in  the  sense  that  the  amount  of  gap  decrease  will  be  smaller  at  each 
iteration  (as  a  consequence  of  the  increasing  values  of  the  centering  parameters  cr's). 

Once  condition  (9)  is  satisfied,  the  value  of  p  will  be  decreased  and  the  gap  will  decrease 
by  a  factor  close  to  or  equal  to  the  value  of  <r°,  in  two  contiguous  iterations. 

The  process  will  continue  in  this  way  with  the  interesting  feature  that  the  satisfaction  of 
condition  (9)  will  force  the  iterates  to  stay  close  to  the  central  path  when  the  gap  is  small. 
Then,  condition  (7)  from  Lemma  4.1  will  be  satisfied  for  a  subsequence  of  the  iteration 
sequence  {(xk,yk,zk)}  generated  by  the  algorithm. 

6  Numerical  Experience 

In  this  section  we  discuss  the  numerical  results  obtained  by  using  the  Long-Step  Shrinking- 
Neighborhood  Algorithm  for  DEA  problems. 

The  experiments  were  performed  in  64  bit  arithmetic  using  a  code  implemented  in 
MATLAB1.  The  code  is  implemented  following  Lustig,  Marsden  and  Shanno  [10],  so  the 
starting  point  is  not  necessarily  feasible.  The  code  generates  a  sequence  of  iterates  that 
approach  feasibility  and  drive  the  gap  to  zero. 

1  Trademark  of  The  Math  Works,  Inc.,  21  Eliot  Street  South  Natick,  MA  01760. 
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We  say  that  a  problem  is  solved  to  an  accuracy  of  10  m  for  some  positive  integer  m  if 
the  algorithm  is  terminated  when 


(\ctxk-bTyk\  \\Axk  -  6||!  \\ATyk  +  zk  -  c\U  pT*z*  -  e||2^ 

V  !  +  l&V|  ’  1  +  Nli  ’l  +  Mi  +  ||^||i’  xTz/n 


In  this  study  all  the  problems  were  solved  to  an  accuracy  of  10  8  using  a  Sun  4/490  work¬ 
station.  The  algorithm  stops  when  the  problem  is  solved  to  the  given  accuracy  or  when  the 
number  of  iterations  is  equal  to  200. 


6.1  Classification  of  the  D Mil's 

There  are  two  main  objectives  in  this  section.  The  first  one  is  to  show  that  the  solution 
computed  by  using  the  Long-Step  Shrinking-Neighborhood  Algorithm  is  the  analytic  center 
solution  for  each  of  the  considered  problems.  The  second  objective  is  to  use  the  computed 
solution  and  the  following  Lemma  (Lemma  10 A  from  [4])  in  order  to  classify  the  DMU's  of 
each  problem  into  the  six  classes  E,  E',  F,  NE,NE',NF. 

Lemma  6.2  (Charnes,  Cooper,  Thrall)  Let  (0O,  X 0,S0)  and  (uo,v0,to)  be  a  SCSC  solu¬ 
tion  of  the  problems  (1)  and  (2),  for  a  selected  DMU0. 

Then,  DMUa  belongs  to 

E  iff  0O  =  1,  S0  =  0,  and  toj  >  0  for  all  j  /  o. 

E'  iff  6 o  =  1,  S0  =  0,  and  tQj  —  0  for  some  j  ^  o. 

F  iff0o  =  l,So=£0. 

N E  iff  0o  <  1,  SQ  =  0,  and  for  some  DMUk  XQk  7^  0  and  A aj  =  0  for  all  j  ^  k. 

N E  iff  60  <  1,  S0  =  0,  and  more  than  one  A 0j  is  nonzero. 

NF  iff  0o  <  l,  S0  =£  0. 
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The  problems  considered  come  from  Thompson,  Langeimeier,  C.  Lee,  E.  Lee  and  Thrall 
[17]  and  Thompson,  Lee  and  Thrall  [18]. 

We  present  a  summary  of  the  results  obtained  for  the  following  problems: 

(PI)  Kansas  Farming  Problem:  Problem  I  from  [17].  This  problem  consists  of  32  wheat 
farms  (32  DMU's),  4  inputs  and  1  output,  with  normalized  output. 

(P2)  Oil-Gas  Firms  Problem:  Problem  II  (Ratio  case)  from  [18].  This  problem  consists  of 
30  Oil-Gas  Extraction  firms  (30  DMU's),  4  inputs  and  2  outputs. 

Following  the  notation  used  in  the  paper,  we  denote  by  ( 9 ,  A,  S )  and  (u,v,t)  solutions  for 
the  linear  programming  problems  (1)  and  (2)  respectively. 

As  an  illustration  of  the  output  obtained  by  using  the  algorithm,  we  show  the  results  for 
DMUi,  DMUS  and  DMU32  from  problem  PI,  in  Tables  1,  2  and  3.  By  counting  the  number 
of  nonzero  values  it  may  be  observed  that  the  solutions  obtained  are  SCSC  solutions  in  each 
problem,  so  any  of  these  Tables  may  be  directly  used  to  classify  the  respective  DMU's  using 
Lemma  6.2. 

Tables  4,  7  contain  the  summarized  information  needed  for  Lemma  6.2  in  order  to  classify 
the  different  DMU's  for  problems  PI  and  P2  respectively. 

Tables  5,  8  contain  the  dual  or  multiplier  solutions  for  the  DMU's  in  E,  for  problems 
PI  and  P2  respectively. 

Tables  6,  9  contain  the  nonzero  values  of  the  primal  or  envelopment  solutions  for  all  the 
DMU's  for  problems  PI  and  P2  respectively. 

Columns  2  and  3  of  Tables  4,  7  show  that  the  solutions  obtained  for  each  of  the  DMU's 
are  SCSC  solutions. 

For  the  Kansas  Farming  Problem  (problem  PI)  we  obtain: 

E  =  { DMUj ,  forj  =8,9,14,15,16,31} 

NE'  =  {DMU32} 
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All  other  DMU's  belong  to  NF. 


For  the  Oil-gas  Firms  Problem  (problem  P2)  we  obtain: 

E  =  {DMUj,  for  j  =  1,2,3,5,9,11,12,15,19,23,24,26,29} 

All  other  DMU's  belong  to  NF. 

6.2  Sensitivity  Analysis 

For  each  extreme  efficient  DMU0  to  find  a  SCSC  solution  centrally  positioned  in  the  mul¬ 
tiplier  space  W0  is  of  extreme  importance  in  evaluating  data  irregularities  (see  Thompson, 
Dharmapala,  and  Thrall  [16]  for  a  detailed  discussion)  . 

The  Long-Step  Shrinking- Neighborhood  Algorithm  presented  in  this  work  was  designed 
for  computing  a  SCSC  solution  being  a  center-most  solution  in  the  sense  defined  in  sec¬ 
tion  3.  This  is  the  first  time  that  a  centrality  condition  is  considered  for  solving  DEA  linear 
programming  problems. 

Thompson,  Dharmapala,  Diaz,  and  Thrall  [15]  studied  the  effect  of  using  the  solutions  ob¬ 
tained  by  the  Long-Step  Shrinking-Neighborhood  Algorithm  to  analyze  the  data.  They  made 
a  multiplier  sensitivity  analysis  to  the  data  for  the  Oil-Gas  Firms  Problem  (Problem  (P2)) 
considering  the  analytic  center  solutions  from  Table  7  and  the  SCSC  solutions  calculated  us¬ 
ing  the  simplex  algorithm  repeatedly  (simplex  derived  centers,  see  [4]).  They  observed  that 
the  solutions  obtained  by  using  the  Long-Step  Shrinking-Neighborhood  Algorithm  improved 
the  stability  of  the  sensitivity  results.  This  study  suggests  that  the  centrality  condition  sat¬ 
isfied  by  the  analytic  center  solutions  can  be  used  to  contribute  additional  robustness  in  the 
DEA  theory. 
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Table  1:  RESULTS  FOR  THE  KANSAS  FARM  PROBLEM:  DMU1 


s 

e 

A 

w 

VTXo  —  1 

t 

0.477352519e-ll 

0.731225593611 

0.527656801e-03 

0.304820975e-10 

0.119477025e-ll 

0.63486621  le-f  00 

0.829983403e-ll 

0.829983403e-ll 

0.133743123e-ll 

0.5024912636-10 

0.469009492e-ll 

0.849414685e-12 

0.764230363e-ll 

0.789475664e-ll 

0.113349502e+00 

0.419171978e-10 

0.441933086e-ll 

0.579963577^11 

0.381 167080e-ll 

0.926337792e-ll 

0.2576681826-10 

0.110093286e+00 

0. 77655721 2e+00 

0.149825333e-ll 

0.500714380e-ll 

0.4388557486-11 

0.190747980e-ll 

0.318297036e-ll 

0.191582653e-ll 

0.274818353e-ll 

0.417850005e-ll 

0.196422959e-10 

0.927996340e-12 

0.425721304e-ll 

0.515330932e-ll 

0.526099609e-ll 

0.127342959e-ll 

0.552632292e-ll 

0.821066350e-ll 

0. 63 486621 le-(- 00 

0.414448000e+00 

0.574341096e-08 

0.994206469601 

0. 253651 265e-|-01 

0.477352519611 

0.365133789e+00 

0.365133789e+00 

0.22659481 7e+01 

0.603104984601 

0.6461 59598e+00 

0.356780958e+01 

0.396549260e+00 

0.383868685e+00 

0.267363314610 

0.722984839e-01 

0.685748576e+00 

0.522541409e+00 

0.795071 246e+00 

0.3271 53861e+00 

0.117614438e+00 

0.275271086e-10 

0.390254550e-ll 

0.202272191e+01 

0.605245220e+00 

0.6905571 74e+00 

0.158877166e+01 

0.952 113750e+00 

0.158184982e+01 

0.110274653e+01 

0.725272183e+00 

0.154286946e+00 

0.326569160e+01 

0.711862390e+00 

0.5880783906+00 

0.576041077e+00 

0.237983308e+01 

0.548384504e+00 

0.369099263e+00 

22 


Table  2:  RESULTS  FOR  THE  KANSAS  FARM  PROBLEM:  DMU8 


s 

6 

A 

w 

VTXo  —  1 

t 

0.570105636D-13 

0.89797X989D-13 

0.459608832D-13 

0.132356100D-12 

0.24605 7010D-1 2 

0.100000000D+01 

0.440773806D-13 

0.753690564D-14 

0.86941871 5D-13 

0.299216371D-13 

0.896999832D-14 

0.330838528D-13 

0.6531 23867D-13 

0.100000000D+01 

0.85361 2776D-13 

0.303352400D-13 

0.446189993D-13 

0.343073471D-13 

0.337956822D-13 

0.706490508D-13 

0.459792663D-12 

0.100228656D-12 

0.940300020D-14 

0. 536671 569D- 13 

0.624274326D-13 

0.133021584D-13 

0.341685228D-13 

0.111535134D-13 

0.319042607D-13 

0.22590041  ID-13 

0.798952502D-13 

0.978546423D-14 

0.415468840D-13 

0.267488973D-13 

0.281631954D-13 

0.684349345D-1 4 

0.205195777D-13 

0.4491 15 117D- 13 

0.100000000D+01 

0.634881314D+00 

0.124041489D+01 

0.430736201D+00 

0.231696564D+00 

0.570105636D-13 

0.129341995D+01 

0.756418700D+01 

0.655731958D+00 

0.190532902D+01 

0.635569391D+01 

0.172321416D+01 

0.872890527D+00 

0.570105636D-13 

0.667873833D+00 

0.187935100D+01 

0.127771946D+01 

0.166175961D+01 

0.168691856D+01 

0.806954417D+00 

0.1 23991 895D+00 

0.568805030D+00 

0. 606301 845D+01 

0.106229894D+01 

0.913229349D+00 

0.428581301D+01 

0. 166851 122D+01 

0.511144440D+01 

0.178692633D+01 

0.252370341D+01 

0.713566370D+00 

0.582604589D+01 

0.137219830D+01 

0.213132388D+01 

0.202429316D+01 

0.833062296D+01 

0.277834976D+01 

0.126939757D+01 

23 


Table  3:  RESULTS  FOR  THE  KANSAS  FARM  PROBLEM:  DMU32 


s 

0 

A 

w 

< 

H 

H 

0 

1 

1 

0.285312080D-14 

0.734008859D-14 

0.258825573D-15 

0.273494886D-13 

0.125575227D-14 

0.647165747D+00 

0.488299879D-14 

0.11058994 7D-14 

0.156970572D-13 

0.294661638D-14 

0.486222149D-15 

0.309173120D-14 

0.696961 199D-14 

0.451224452D-01 

0.158130335D-13 

0.310002443D-14 

0.421937525D-14 

0.340843431D-14 

0.793726566D-14 

0.126394816D+00 

0.235073169D+00 

0.593409570D+00 

0.127445765D-14 

0.384180336D-14 

0.362425758D-14 

0.159914870D-14 

0.237814921D-14 

0.128104662D-14 

0.166311639D-14 

0.310344576D-14 

0.808334773D-14 

0.701951308D-15 

0.348980844D-14 

0.728381201D-15 

0.266476701D-14 

0.108433758D-14 

0.478037926D-14 

0.523317122D-14 

0.647165747D+00 

0.251555827D+00 

0.713392435D+01 

0.675128551D-01 

0.147038720D+01 

0.28531 2080D-14 

0.378136907D+00 

0.166962921D+01 

0.117629823D+00 

0.626631302D+00 

0.379752765D+01 

0.597219465D+00 

0.264927525D+00 

0.409207003D-13 

0.116767099D+00 

0.595621777D+00 

0.437610297D+00 

0.541727341D+00 

0.232629489D+00 

0.146085268D-13 

0.785475462D-14 

0.311158119D-14 

0.144880613D+01 

0.48061 8575D+00 

0.509467667D+00 

0.115464062D+01 

0.776419767D+00 

0.144135430D+01 

0.1 1 1023021D+01 

0.594965145D+00 

0.228425415D+00 

0.2630441 79D+01 

0.529095534D+00 

0.25349941 1D+01 

0 .692909381 D +00 

0.170282953D+01 

0.386254302D+00 

0.352834253D+00 

24 


Table  4:  RESULTS  FOR  THE  KANSAS  FARM  PROBLEM 


DMU 

6 

mixij  |a,  +  tojj 

mi 

{}  ■  Aj  >  0} 

CLASS 

1 

6.3487e-01 

5.2766e-04 

6.0310e-02 

5.2766e-04 

8  15  16 

NF 

2 

4.5162e-01 

2.3815e-02 

2.1789e-02 

5.4174e+00 

8  14 

NF 

3 

9.6969e-01 

4.3558e-03 

3.0309e-02 

4.3558e-03 

9  15  16 

NF 

4 

5.3024e-01 

3.4651e-01 

2.3649e-02 

4.1843e-01 

8  14  16 

NF 

5 

1.5106e-01 

3.4519e-03 

8.4651e-03 

3.4519e-03 

8  15  16 

NF 

6 

7.8899e-01 

2.9878e-02 

4.7398e-02 

1.8294e+00 

8  16 

NF 

7 

7.7507e-01 

3.3121e-02 

3.7395e-02 

1.9713e-01 

8  14 

NF 

8 

l.OOOOe+OO 

2.3170e-01 

1.2399e-01 

5.7118e-13 

8 

E 

9 

1.0000e+00 

2.6241e-02 

7.9946e-02 

4.6025e-12 

9 

E 

10 

8.1561e-01 

1.1183e-02 

1 ,7065e-01 

3.1344e+00 

8 

NF 

11 

6.3743e-01 

7.6213e-03 

3.0754e-02 

4.9780e-01 

8  14 

NF 

12 

6.8896e-01 

2.7369e-02 

5.041 7e-02 

5.5087e-01 

14  15 

NF 

13 

9.4779e-01 

6.1937e-01 

1.6970e-02 

2.7653e+00 

8  14  16 

NF 

14 

l.OOOOe+OO 

1.7030e-02 

2.0229e-01 

4.4938e-12 

14 

E 

15 

l.OOOOe+OO 

2.8111e-02 

6.0560e-01 

1.2408e-12 

15 

E 

16 

— 

l.OOOOe+OO 

2.9851e-02 

7.2386e-02 

4.7218e-12 

16 

E 

17 

4.2907e-01 

2.2462e-02 

1.8105e-01 

2.7632e+00 

14 

NF 

18 

6.431  le-01 

3.9577e-02 

4.7062e-02 

8.5941e-02 

14  15 

NF 

19 

6.8578e-01 

5.1246e-02 

5.0184e-02 

1.4179e-01 

14  15 

NF 

20 

4.3438e-01 

4.0239e-02 

8.9820e-03 

1.1709e+00 

8  14 

NF 

21 

5.0561e-01 

4.9383e-02 

7.3112e-02 

4.9383e-02 

8  14  15 

NF 

22 

3.6053e-01 

2.3561e-01 

1.9861e-02 

1.5674e+00 

8  14  16 

NF 

23 

3.7680e-01 

4.9142e-04 

2.6901e-01 

8.3244e-03 

8  15 

NF 

24 

6.0339e-01 

3.9431e-01 

3.3238e-02 

1.5156e+00 

8  14  16 

NF 

25 

8.0449e-01 

8.0235e-03 

5.4330e-02 

8.0235e-03 

8  15  16 

NF 

26 

2.3553e-01 

3.6162e-02 

1.7236e-02 

1.6238e-01 

14  15 

NF 

27 

6.8683e-01 

— 

1.1417e-02 

5.0261e-02 

5.8361e-01 

14  15 

NF 

28 

7.4989e-01 

3.9878e-02 

3.6180e-02 

7.8296e-01 

8  14 

NF 

29 

5.5820e-01 

1.0230e-02 

3.3534e-02 

4.2322e-01 

8  16 

NF 

30 

4.5968e-01 

2.1411e-02 

2.2178e-02 

5.9699e+00 

8  14 

NF 

31 

l.OOOOe+OO 

9.4731e-05 

9.5463e-04 

1.8868e-ll 

31 

E 

32 

6.471 7e-01 

6.7513e-02 

4.5122e-02 

3.9057e-14 

8  14  15  16 

NE’ 

25 


Table  5:  RESULTS  FOR  THE  KANSAS  FARM  PROBLEM 


DUAL  MULTIPLIERS  w  =  (u,  «)  FOR  THE  DMU'i  IN  E 

DMU 

U1 

VI 

V2 

V3 

V4 

8 

1.0000e  +  00 

6.3488e-01 

1.2404e-f  00 

4.3074e-01 

2.31 70e-01 

9 

1.0000e  +  00 

3.4861e-02 

5.9£72e-01 

2.6241e-02 

2.8584e+01 

14 

1.0000e+00 

3.9452e-02 

2.099Se+01 

1.7030c- 02 

1.1569e+01 

15 

1.0000e+00 

2.8111e-02 

1. 0096c -f  00 

1.0009e+00 

2.5996e+00 

16 

1.0000e+00 

4.8521e-01 

2.2580e-f  00 

2.9851e-02 

1.4103e+01 

31 

1.0000e+00 

1.9060e-01 

4.1392e+00 

9.4731e-05 

2.0512e+01 

Table  6:  RESULTS  FOR  THE  KANSAS  FARM  PROBLEM 


DMU 

NONZERO  PRJMAL  VARIABLES  X 

1 

Xg  =  1 ,1335c  -  01 

Xj5  =  1.1009c  -  01 

A16  =  7.7656c  -  01 

_ 

2 

- 

- 

3 

Ag  =  7.3275c  -01 

A14  sb  1. 7646c  -  01 

Xjg  =  9.0795c  -02 

- 

4 

A8  =  3.7650c  -01 

Xig  =  5.9985c  -01 

- 

5 

Xj4  =  8.2373c  -01 

X16  =  8.4651c  -  03 

- 

6 

Ag  =  4.0857e  -  01 

X16  =  5.9143c  -01 

- 

_ 

7 

Xg  =  S. 3739c  -01 

Xj4  =  4.6261c  -01 

- 

- 

8 

Ag  =  1.0000c  +  00 

- 

- 

- 

9 

Xg  =  1  .OOOOe  +  00 

- 

- 

- 

10 

Xg  =  1. OOOOe  +  00 

- 

- 

- 

11 

Ag  =  7.0385e  -  01 

X14  =  2. 9615c  -  01 

- 

- 

12 

x14  =  5. 4344  c  -  01 

Xls  =  4. 5656c  -  01 

- 

- 

13 

Ag  =  3.8210c  -  01 

Xj4  =  1.6970c  -  02 

- 

14 

A14  =  1  .OOOOe  +  00 

- 

- 

- 

15 

- 

- 

- 

16 

X16  =  1.0000c  +  00 

- 

- 

- 

17 

X14  =  1.0000c  +  00 

i 

- 

- 

18 

Xj6  =  7.5631c  -  01 

- 

- 

19 

- 

- 

20 

Aj  =  4.3438c  -  01 

Xg  =  8.9820c  -  03 

Xu  =  9.9102c  -  01 

- 

21 

Xg  =  2.5511c  -  01 

Xj4  =  3.5935c  —  01 

Xj5  =  3.8554c  -  01 

- 

22 

Xg  =  1.0156c  —  01 

A14  =  7.2329c  -02 

- 

23 

Ag  =  6.1535c  -  01 

Xj5  =  3. 8465c  -  01 

- 

24 

Xg  =  2.9293c  -  01 

- 

25 

Ag  =  5.4330c  -  02 

A15  =  3.8062c  -  01 

Aj6  =  5. 6505c  -  01 

- 

26 

Xj4  =  5. 9808c  -  01 

KSBESESBI 

- 

- 

27 

X14  =  2.0409c  -  01 

X15  =  7.9S91e  -01 

- 

- 

28 

Xg  =  8.8701c  —  02 

X14  =  9.1 130e  -  01 

- 

- 

29 

Xg  =  1.2523c  —  01 

A16  =  8.7477c  -  01 

- 

_ 

30 

Ag  =  1.0537e  -  01 

Xj4  =  8.9463c  -  01 

- 

- 

31 

X31  =  1  .OOOOe  +  00 

- 

- 

- 

32 

Xj5  =  2.3507c  -  01 

A16  =  5.9341e  -01 

26 


Table  7:  RESULTS  FOR  THE  OIL/GAS  FIRMS  PROBLEM 


DMU 

e 

minj  |sj  +  wj\ 

KfflfflVSH 

lUlli 

{j  :  Xj  >  0} 

CLASS 

1 

1.0000e+00 

1.1395e-04 

1.4141e-02 

7.5891e-ll 

1 

E 

2 

l.OOOOe+OO 

2.1788e-03 

4.5875e-02 

5.4520e-13 

2 

E 

3 

l.OOOOe+OO 

2.5435e-03 

3.3921e-03 

4.6569e-ll 

3 

E 

4 

7.2614e-01 

1.1798e-01 

7.8232e-03 

1.5079e+00 

9  23  26 

NF 

5 

l.OOOOe+OO 

3.0406e-03 

2.3723e-01 

2.3705e-14 

5 

E 

6 

9.4024e-01 

1.4180e-03 

3.9281e-03 

5.9584e+00 

9  19  26 

NF 

7 

7.9612e-01 

4.1168e-04 

1.2080e-02 

3.7405e+00 

15  19  26  29 

NF 

8 

5.5799e-01 

7.8339e-04 

4.3760e-02 

4.4637e-01 

11  19  23  26 

NF 

9 

1.0000e+00 

1.1820e-03 

1 ,5607e-02 

1.2403e-ll 

9 

E 

10 

6.4746e-01 

6.5483e-03 

2.5451e-03 

7.1769e-02 

1  5  24  26 

NF 

11 

l.OOOOe+OO 

1.7149e-04 

4.4009e-04 

8.2657e-18 

11 

E 

12 

1.0000e+00 

2.4588e-03 

2.0083e-02 

7.4641e-15 

12 

E 

13 

6.0141e-01 

6.7137e-04 

1.0454e-02 

1.2281e+00 

15  19  26  29 

NF 

14 

2.6143e-01 

8.7400e-03 

2.7892e-03 

1 .0903e+00 

1  26 

NF 

15 

l.OOOOe+OO 

1.6835e-03 

2.5313e-01 

5.7282e-14 

15 

E 

16 

3.9684e-01 

8.6391e-03 

4.2020e-03 

9.4331e-01 

11  19  23  26 

NF 

17 

7.9127e-01 

3.2060e-03 

6.6778e-03 

5.3531e-01 

9  19  26 

NF 

18 

6.1208e-01 

1.5648e-02 

4.7808e-03 

1.5648e-02 

2  19  23  26 

NF 

19 

l.OOOOe+OO 

4.4342e-02 

4.6668e-01 

6.1287e-15 

19 

E 

20 

5.4550e-01 

2.3840e-02 

1.5682e-02 

7.4206e-02 

2  19  23  26 

NF 

21 

5.7068e-01 

9.0396e-03 

6.3437e-03 

4.3583e-01 

9  26 

NF 

22 

9.1175e-01 

1.3689e-02 

1.1081e-02 

3.0801e+00 

15  26  29 

NF 

23 

l.OOOOe+OO 

2.7210e-02 

l.OOOOe+OO 

1.0015e-15 

23 

E 

24 

l.OOOOe+OO 

2.6999e-05 

5.3297e-03 

4.0545e-14 

24 

E 

25 

8.5307e-01 

1.6499e-03 

4.9890e-03 

3.3016e+00 

2  19  23  26 

NF 

26 

l.OOOOe+OO 

1.6248e-03 

3.2859e-02 

1.7429e-14 

26 

E 

27 

4.6902e-01 

1.8874e-02 

1.1944e-02 

1.8138e+00 

15  26  29 

NF 

28 

5.3820e-01 

3.8300e-03 

8.3330e-03 

9.1593e-02 

11  19  23  26 

NF 

29 

l.OOOOe+OO 

1.1265e-03 

2.6538e-02 

1.4605e-14 

29 

E 

30 

6.7486e-01 

4.4651e-03 

4.4619e-03 

4.7914e-01 

15  26  29 

NF 

27 


Table  8:  RESULTS  FOR  THE  OIL/GAS  FIRMS  PROBLEM 


DUAL  MULTIPLIERS  ti>  as  (u,  v)  FOR  THE  DMU’«  IN  B 

DMU 

U1 

U2 

VI 

V2 

V3 

V4 

1 

1.0053e-01 

7.5247e-02 

1.1407e-02 

8.4621e-02 

1 . 1395e-04 

5.8077e-02 

2 

5.2761e-01 

3.59S3e-01 

1.6787e-02 

1.3186e+00 

2. 17880-03 

7.0068c- 02 

3 

2.l734e+00 

3.7393e-02 

3.3691e-03 

2.9802O-01 

3.3934e-02 

2.5435e-03 

4 

1.2225e  +  00 

2.6358e-01 

4. 7415o-03 

3.0406O-03 

3. 77620-03 

8.3022e-01 

5 

9.8530e+00 

7.1545e-02 

2.0850e-03 

1.947«o+00 

7.91 75o-02 

1 . 1820e-03 

6 

7.6593e+00 

1.7115e+00 

6. 43740-01 

1.1002c +00 

2.4027e-01 

1 . 7149e-04 

7 

2.6326e+00 

2.4588e-03 

2.5770O-03 

7.0159e-03 

3.8157e-02 

2.7267e-01 

8 

3.8882e  +  00 

3.1 194e-01 

7.4720O-03 

2.7055e-02 

9.6971e-01 

1.6835e-03 

9 

7.7356e+00 

2.5336e+00 

4. 4342o-02 

3.2403e+00 

1. 5156o-01 

1.0898o+00 

10 

6.9897e+01 

8.9983e+00 

3.7751e-01 

2.B099O+02 

5.8142e-02 

2.7210e-02 

11 

3.1337e-01 

5.5737e-04 

3.5389e-03 

4.1623c-04 

2.6999e-05 

7.3642e-02 

12 

7.6442e-01 

2.8138e-02 

1.7083e-01 

1.0578c-02 

1.6248e-03 

3.4998e-02 

13 

1.1507e+01 

2.5467e-01 

4.5919e-02 

1.1265e-03 

1.1704O+00 

9.3267e-02 

Table  9:  RESULTS  FOR  THE  OIL/GAS  FIRMS  PROBLEM 


DMU 

NONZERO  PRIMAL  VARIABLES  A  | 

1 

Aj  =  1 ,0000c  +  00 

- 

- 

- 

2 

A2  sb  1.0000c  +  00 

- 

- 

- 

3 

A3  ss  1.0000c  +  00 

- 

- 

- 

4 

A9  =  7.8232e  -  03 

A23  =  2.6942e  +  00 

A26  =  2.6125c  -  02 

5 

A5  =  1.0000c  +  00 

- 

- 

- 

6 

Ag  =  5.4611c  -  01 

A19  =  4.8105c  -  01 

A26  =  1.1783e  -  01 

- 

7 

A15  =  8.9817e  -02 

A19  =  1.1871«  +  00 

A26  as  1.4027c  -01 

A29  =  2.5732c  -02 

8 

Jkj!  -  1.6993c  -  01 

Ajg  =  1.0968c  +  00 

A23  =  8.3803c  -  01 

A26  =  4.7418c  -  02 

9 

Ag  =  1 .0000  c  +  00 

- 

- 

- 

10 

Aj  =  2.5451c  -  03 

A5  =  1.1024c  -  01 

A24  =  5.1342c  -  03 

A26  =  2.4428c  -  02 

11 

Ajj  =  l.OOOOe  +  00 

- 

- 

- 

12 

A12  =  l.OOOOe  +  00 

- 

- 

- 

13 

A19  =  7.9042c  -  01 

A26  =  5.1826c  -  02 

A29  as  1.0454c  -  02 

14 

Aj  =  2. 7892c  -  03 

15 

Aj  =  l.OOOOe  +  00 

- 

- 

- 

16 

An  =  1.9063c  -  02 

A19  =  4.6213c  -  02 

A23  =  3.0811c  -  02 

A26  =  4.2020c  -  03 

17 

Ag  as  9.4548c  -  03 

A26  =  6.6778c  -  03 

- 

18 

A2  =  4.7808c  —  03 

Ajg  =  8.0622e  —  02 

A23  as  2.0079c  -  01 

A26  =  9.2983c  -  03 

19 

A19  =  l.OOOOe  +  00 

- 

- 

- 

20 

A2  =  1.2193e  -  01 

Ajg  =  1.7270c  -  01 

A23  =  5.1396c  -  01 

A26  =  1.5682e  -  02 

21 

A9  =  1.4964c  -  01 

A26  =  6.3437c  —  03 

- 

- 

22 

A15  as  4.0172e  -  02 

A26  =  1.1081e  -  02 

A29  =  1.7340c  -  0223 

A23  as  l.OOOOe  +  00 

24 

A24  as  l.OOOOe  +  00 

- 

- 

- 

25 

A2  =  1.3020c  +  00 

A19  =  9.4998e  +  00 

A23  =  2. 3073c  +  01 

A26  =  5.5866c  —  01 

26 

A26  =  1  OOOOe  +  00 

- 

- 

- 

27 

A15  =  2.3529e  -02 

A26  =  1.1944c  -  02 

A29  =  4. 8987c  -  02 

- 

28 

mirmimmiM 

A19  =  6.6127c  -  02 

A23  —  7.3328e  -  01 

A2(j  =  8.3330c  —  03 

29 

A2g  =  l.OOOOe  +  00 

- 

- 

- 

30 

A26  =  1  0155c  -  02 

A29  =  4.4619e  -  03 

- 

28 


T  Concluding  Remarks 


In  this  paper  a  primal-dual  interior-point  algorithm  for  solving  the  DEA  linear  programming 
problems  was  proposed.  The  algorithm  is  a  modification  of  the  Kojima-Mizuno-Yoshise 
primal-dual  algorithm.  The  sequence  of  iterates  generated  by  the  algorithm  converges  to  the 
solution  set  while  approaching  the  central  path;  hence  the  computed  solution  is  the  analytic 
center  of  the  solution  set.  The  approach  to  the  central  path  is  done  in  such  a  way  that  long 
steps  can  be  taken  when  the  gap  is  not  small,  but  close  to  the  solution  set  only  small  steps 
are  allowed  because  of  the  use  of  gradually  shrinking  neighborhoods  of  the  central  path. 

Numerical  results  were  presented  for  some  real  world  DEA  problems  in  order  to  illustrate 
the  performance  of  the  algorithm.  The  solution  computed  by  the  algorithm  can  be  effectively 
used  for  the  classification  of  the  DMU's  for  DEA  problems. 

The  effective  use  of  the  LemmalOA  from  [4]  for  classifying  the  DMU's  strongly  depends 
on  the  computation  of  a  SCSC  solution  being  far  of  the  relative  boundary  of  the  solution  set, 
in  the  sense  introduced  in  Section  2.  Since  the  solution  computed  by  the  algorithm  proposed 
in  the  present  work  is  a  center-most  solution  in  the  sense  that  it  maximizes  the  product  of 
the  positive  components  among  all  the  SCSC  solutions,  it  seems  to  be  a  good  choice  for  the 
DEA  purposes,  although  it  may  not  be  the  best  solution  choice.  To  our  knowledge,  this  is 
the  first  time  that  a  centrality  requirement  is  considered  for  finding  solutions  of  the  DEA 
linear  programming  problems.  Further  research  is  needed  to  understand  the  relationship 
between  the  solution  here  obtained  and  the  DEA-center  solution  introduced  in  Section  2, 
as  well  as  the  relationship  between  the  optimal  value  dc0  and  d*  =  d0(w*)  with  w*  the  dual 
multiplier  of  the  analytic  center  solution  and  da  the  sensitivity  function  also  introduced  in 
Section  2. 
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